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UNSTEADY  TURBULENT  BOUNDARY  LAYERS 
IN  ADVERSE  PRESSURE  GRADIENTS 


Eugene  E.  Covert* 

Peter  F.  Lorber** 

Department  of  Aeronautics  and  Astronautics 
Massachusetts  Institute  of  Technology 
Cambridge,  Massachusetts  02139 


Abstract 


A  number  of  characteristics  of  an  unsteady  tur¬ 
bulent  boundary  layer  have  been  measuredon  the 
after  part  of  the  upper  surface  of  an  0012  airfoil 
(x/c  "  .69  and  .94).  The  data  is  taken  at  chord 
Reynolds  number  of  700,000  and  over  a  range  of  re¬ 
duced  frequencies  (based  upon  the  semichord)  of  0.5 
to  6.4.  Mean  and  unsteady  velocity  profiles,  as 
well  as  Reynolds  stress  profiles,  are  presented,  as 
is  the  non-linear  coupling  from  the  unsteady  motion 
into  the  steady  motion.  Data  is  presented  and  dis¬ 
cussed  which  shows  that  for  this  experiment  the 
periodic  unsteady  turbulent  velocity  profile  tends 
toward  a  universal  shape,  independent  of  the  mean 
adverse  velocity  gradient,  at  higher  reduced 
frequencies. 


Nomenclature 


a 

b 

c 


V 

x 


6 

r 

s' 

9 

V 

p 


elliptic  cylinder  minor  axis 
elliptic  cylinder  major  axis 
airfoil  chord 


skin  friction  coefficient  t,./4  P  U  2 

/.  00 

pressure  coefficient,  (P-P  P  U  2 

00  2  *°° 

boundary  layer  shape  parameter,  S  /0 

reduced  frequency,  <oc/2U 

pressure 

Reynolds  number  based  on  momentum  thickness 

u  e/v 


external  velocity 


free stream  velocity 


velocity  component  parallel  to  surface 
friction  velocity,  /td/q 

velocity  component  normal  to  surface 
coordinate  parallel  to  surface 
coordinate  normal  to  surface 
airfoil  angle  of  attack 

pressure  gradient,  6/ 

dx 

pressure  gradient,  ~~  (Re„)1/4 

ue  aX  17 

oo 

displacement  thickness, /  (1-u/U  )  dy 

J  o  e 

oo 

momentum  thickness,  f  u/U  (1-u/U  )  dv 
/  e  e  ' 

J  o 

kinematic  viscosity 
density 

shear  stress  at  the  surface 
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6)  radian  frequency 

__  time  average 

'v  period  component  _  % 

<  >  ensemble  average,  <u>  =  u  +  u 

■  nonperiodic  fluctuating  component  of  zero 

mean,  u  =  <u>  +  u* 

Introduction 

Recently  L.  W.  Carr^presented  a  survey  of  the 
state-of-the-art  of  measurements  on  unsteady  tur¬ 
bulent  boundary  layers.  In  many  of  these  experi¬ 
ments  the  metric  surface  is  a  flat  plate  of  one  wall 
of  the  wind  tunnel.  Unsteadiness  is  frequently 
introduced  by  oscillating  vanes  or  shutters  located 
either  upstream  or  downstream  of  the  test  surface. 

Of  all  the  experiments  cited,  those  of  Karlsson  , 
Patel ' ,  Kenison  ,  Schachenmang  and  Rockwell  , 
Simpson,  Shivaprasad  and  Chew  ,  and  Cousteix, 
Houdeville  and  Raynaud^  have  produced  data  that  are 
perhaps  most  closely  related  to  that  which  will  be 
presented  below.  Several  years  ago  Telionis®  re¬ 
viewed  the  state  of  understanding  of  both  separated 
and  attached  unsteady  boundary  layers.  He  concluded 
that, 

a)  generally  the  mean  profiles  are  affected 
very  little  by  the  unsteadiness, 

b)  the  pressure  gradient  has  a  strong  effect 
on  the  fluctuations,  and 

c)  the  state-of-the-art  of  turbulence  model¬ 
ling  is  inadequate. 

An  example  of  the  latter  is  calculations  that  pre¬ 
dict  an  overshoot  in  unsteady  velocity  amplitude  of 
1-2%  over  that  at  the  edge  of  the  boundary  layer, 
while  the  data  shows  overshoots  of  10  times  that 
value  for  the  same  conditions.  Telionis  concluded 
something  fundamental  was  missing  from  the  models. 

In  this  paper,  we  will  present  additional  data 
that  characterizes  the  unsteady  turbulent  boundary 
layer.  Our  metric  surface  is  the  upper  surface  of 
a  N.A.C.A.  0012  airfoil  which  is  fixed  in  a  wind 
tunnel.  As  described  below  the  unsteady  flc.’  is 
generated  by  aerodynamic  interference  at  the  trail¬ 
ing  edge.  We  will  present  data  on  <u>,  <v>,  <v’^> 
and  <u'v’>  and  offer  some  comments  on  the  processes 
and  their  limits. 


■  otion  of  Experiment 

The  expel1..  re  conducting  is  somewhat 

different  i;  con-  ■  •  i  from  those  refered  to  in  the 
introduction.  As  indicated  in  the  introduction,  our 
metric  surface  is  a  N.A.C.A.  0012  airfoil,  with  chord 
length  50.8  cm.  It  is  located  between  two  vertical 
sidewalls  in  the  MIT  Wright  Brothers'  Wind  Tunnel. 

The  walls  are  not  quite  parallel  to  the  flow,  rather 
they  diverge  enough  to  compensate  for  the  effects  of 
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boundary  layer  growth.  Preliminary  data  taken  by 
Kanevsky9  on  a  flat  plate  with  a  25  micro  inch  sur¬ 
face  installed  in  place  of  the  airfoil  showed  that 
the  CDon  the  plate  was  uniform  to  within  +  1%  over 
the  entire  plate  and  +  .3%  over  the  metric  zone 
(Fig.  1).  The  measured  velocity  profiles  were  ex¬ 
pressed  in  the  form  of  the  law  of  the  wall  as  (Fig.  2) 

u 

—  =  2.5  In  -  +5.5 

U  V 

T 

The  pressure  distribution  on  the  airfoil  sur¬ 
face  due  to  the  unsteady  flow  generated  by  the 
rotating  ellipse  is  described  in  detail  in  Ref.  10 
and  11.  For  convenience  some  results  of  the  pres¬ 
sure  measurements  are  summarized  here.  The  current 
boundary  layer  experiments  were  performed  at  a 
tunnel  velocity  of  20  mps,  resulting  in  a  Reynolds 
number  based  on  airfoild  chord  of  7  x  10^.  At  this 
velocity  the  freestream  turbulence  level  is  approx¬ 
imately  0.7%.  Figures  4  and  5  show  the  mean  pres¬ 
sure  coefficient  distributions  for  this  Reynolds 
number,  for  reduced  frequencies  k  =  0.5  to  6.4, 
and  for  mean  angle  of  attack  of  0  and  10  degrees. 
Figures  6  and  7  illustrate  the  ensemble  averaged 
fundamental  harmonic  pressure  coefficient  amplitude 
and  phase  lag  distributions  seen  at  0  and  10°  angle 
of  attack  for  the  upper  surface  of  the  airfoil. 

Upper  surface  mean  and  unsteady  amplitudes  increase 
smoothly  near  the  trailing  edge,  while  the  phase 
lag  in  this  region  is  nearly  constant.  Note  both 
amplitude  and  phase  distributions  shown  in  Figs.  6 
and  7  were  relatively  independent  of  the  geometric 
angle  of  attack,  as  long  as  the  airfoil  has  not 
stalled. 

12 

Following  Clauser  ,  we  could  define  the  non- 


dimensional  pressure  gradient 


o  dx 
2x  d«e 


,  which  for 


steady  flow  is  equal  to  -  —  ,  a  quantity 

e 

independent  of  Reynolds  number.  The  mean  value  of 
this  parameter  at  the  boundary  layer  measurement 
locations  ranges  for  the  present  test  from  +  0.25 
to  +  1.7.  Clauser  gives  0.58  as  the  maximum  pres¬ 
sure  gradient  for  equilibrium  profiles.  Similarly, 
for  the  unsteady  pressure  distribution,  the  unsteady 

amplitude  of  ^  ranged  from  .02  to  .40.  The 

experimental  conditions  are  given  in  Tables  1  and  2. 

Buri's  form  of  the  mean-nondimensional  pressure 

gradient,  r  =  —  (Re^)1^4  is  also  useful  in 

Ue  dx  9 

defining  test  conditions.  For  our  profiles,  T  was 
between  -.002  and  -.042.  T  a,  -.06  corresponds  to 
separation,  for  a  steady  mean  flow^.  T  may  be 

O  ,  du 

fj  1 

written  in  steady  flow  as  —  which  is 

2  40C  _  U  dx 

f  c 

closely  related  to  Clauser* s  alternate  form  8  = 

6  *  d 

—  ^  .  The  advantage  of  Clauser' s  8  or  Buri's  Y 

w 

is  that  the  boundary  layer  history  is  implicity 
included. * 


Boundary  layer  velocities  were  measured  using 
either  a  single  hot  wire  (determining  longitudinal 
velocity)  or  a  cross  wire  (determining  both  longi¬ 
tudinal  and  normal  velocities).  The  velocity  probes 
were  motor  driven  across  the  boundary  layer,  with 
the  position  being  given  by  a  linear  potentiometer. 
Flow  corporation  constant  temperature  anemometers, 
linear izers  and  sum  difference  amplifiers  were  used 
to  produce  u  and  v  velocities. 

With  the  exception  of  some  of  the  low  pressure 
gj.,_lient  data,  for  which  an  analog  system  was  used, 
the  velocities  and  pressures  were  digitized,  en¬ 
semble  averaged,  converted  to  nondimen sional  form 
and  plotted  on  line  (Fig.  8).  This  procedure 
allowed  monitoring  of  the  progress  of  the  experi¬ 
ment,  and  in  particular  indicated  when  additional 
cycles  of  data  were  needed  in  order  to  properly 
define  the  unsteady  quantities.  The  required  num¬ 
ber  of  cycles  depends on  the  signal  to  noise  ratio 
of  the  periodic  to  fluctuating  amplitudes,  and 
ranges  from  100  to  4096. 

For  da£a  taken  with  the  single  hot  wire,  only 
<u>  =  u  +  u  and  u'u'  were  obtained,  while  data** 
taken  with  the  cross  wire  also  included  <v>,  <u'u'> 
<v'v’>,  <u'v’>,  where  u*  is  defined  to  have  zero 
time  average  and  ensemble  average.  This  data  will 
be  stored  on  tape  and  provided  to  Dr.  L.W.  Carr  at 
Ames  Research  Center1. 


*That  is,  Kantian' s  integral  relation  may  be  used  to 
relate  B  or  Y  to  the  flow  in  the  boundary  layer. 

B  ,  SE  .  JL  (V  30  . 

T  dx  H+2  )u  2  dx  ( 


H_  d0  | 

+2  |ut2  dx  '  ( 

-1  J  V  d0  . 

80 (H+2)  j  Ux2  dx 


Hence  either  Clauser' s  6  or  Buri's  Y  are  measures  of 
a  departures  of  the  turbulent  boundary  layer  from 
the  flat  plate  equilibrium  flow.  If  8  or  f  are  pure 
constants  (ie,  independent  of  Reynolds  number)  then 
the  flow  in  question  is  selfsimilar.  In  any  case  fl 
or  T  represent  the  ratio  of  two  forces;  one  is  force 
per  unit  volume  applied  to  the  boundary  layer  from 
outside  the  boundary  layer  and  the  other  is  the  force 
per  unit  volume  applied  by  the  wall  through  the  shear 
stress,'  where  the  volume  is  normalized  to  either  the 
displacement  or  momentum  defect  thickness.  The 
presence  of  the  products  of  boundary  parameters  im¬ 
plies  the  unsteady  boundary  layer  offers  the  possi¬ 
bility  for  rectification  and  multiples  of  the 
fundamental  frequency  of  excitation. 

**Note  the  decomposition 

_ 

u  =  u  +  u  +  u ' 

only  implies  that  u'  has  a  zero  temporal  mean.  There 
is  no  intent  to  imply  that  u'  is  independent  of  fre¬ 
quency  of  excitation,  ie,  in  the  notation  of  Ref.  14 

u'(x,y,t)  =  (x , y , t )  +  uj  (x,y,t,el,t) 

except  that  we  normalize  such  that  uj  +o  as  t  >o 


TABLE  T.  SUMMARY  01*  EXPERIMENTAL  CONDITIONS 


k 

— 

a 

y/c 

- 

e/c 

=■ 

6 

x  dP 
q  dx 

•v 

X  dP 
q  dx 

mm 

0.5 

ffl  jK 

0.69 

-.0025 

v';,:  H 

1.80 

| 

.301 

.0483 

.021 

1.0 

B9 

1.71 

.281 

.0593 

.022 

1.5 

mm 

.0017 

1.74 

1340 

.422 

.0575 

.014 

2.0 

MM 

.0016 

1.75 

1230 

.334 

.0299 

.0057 

3.9 

.0018 

1.74 

1410 

.297 

.0198 

.0034 

6.4 

I?  ■  | 

-.0021 

.0019 

1.74 

1430 

.255 

.0375 

.0038 

0.5 

2 

0.94 

-.0144 

.0025 

mm 

1.71 

.273 

.031 

mem 

i 

-.0151 

.0026 

1.68 

.402 

.038 

wmm 

-.0138 

.0025 

1.61 

1670 

1.62 

.203 

.013 

mSM 

1 

-.0101 

.0025 

1.61 

1780 

1.19 

.224 

.0078 

M 

10* (1) 

0.94 

-.0408 

.0090 

2.38 

5950 

.969 

.173 

.020 

-.0358 

.0088 

2.50 

5600 

.891 

.228 

.033 

!  2.0 

-.0401 

.0091 

2.18 

5980 

.951 

.000 

.0076 

-.0406 

.0093 

2.04 

6120 

.947 

.089 

.0048 

-.0430 

.0095 

2.00 

6390 

.958 

.109 

.0062 

0.5 

io-(2> 

0.94 

-.0297 

.0093 

2.56 

6390 

^pppplpl 

.238 

.040 

1.0 

-.0322 

.0094 

2.53 

6230 

.262 

.047 

2.0 

-.0297 

.0092 

2.44 

6160 

.689 

.153 

.016 

3.9 

-.0306 

.0087 

2.22 

6420 

.743 

.135 

.0088 

6.4 

_ _ 

-.0304 

.0083 

2.11 

_ _ _ 

6400 

.769 

.166 

.0079 

(1)  Cylinder  axis  at  x/c  »  1.18,  y/c  =  -.28 

(2)  Cylinder  axis  at  x/c  =*  1.12,  y/c  =  -.19 


TABLE  II.  STEADY  STATE  TEST  CONDITIONS 


ELLIPSE 

ORIENTATION 

a 

x/c 

r 

8/c 

% 

H 

x  dP 
q  dx 

HORIZ. 

0” 

.69 

-.0017 

.0013 

970 

1.81 

.55 

VERT. 

-.0032 

.0020 

1500 

1.87 

.57 

HORIZ. 

0»«> 

.94 

-.011 

.0022 

1410 

1.66 

1.58 

VERT. 

-.037 

.0047 

2900 

1.90 

1.67 

HORIZ. 

10-(2> 

.94 

-.034 

.0075 

4910 

2.30 

.80 

VERT. 

-.0078 

.0107 

7300 

3.40 

.13 

U) 


cylinder  axis  at  x/c  »  1.18,  y/c  **  -.28 


(2) 


cylinder  axis  at  x/c  *  1.12,  y/c  »  -.19 


Section  omitted  from  page  2  (after  the  first  equation) 


This  is  in  general  agreement  with  other  flat  plate  turbulent 
boundary  layer  data.  Thus  we  feel  that  the  mean  effects  of 
the  installation  have  been  minimized. 

The  airfoil  may  be  rotated  about  the  trailing  edge  to 
give  a  mean  geometric  angle  of  attack  of  between  0  and  15 
degrees.  Just  downstream  of  and  below  the  the  trailing  edge 
a  two  dimensional  elliptic  cylinder  (b/a  =  2.12)  is  located, 
with  its  axis  at  either  x/c  =  1.175,  y/c  =  -.276  or  x/c  =  1. 
y/c  =  -#19«  This  elliptic  cylinder  is  rotated  at  from  0  to 
3000  r.p.m. ,  through  a  belt  driven  by  an  adjustable  speed 
electric  motor.  This  provides  and  unsteady  perturbation  to 
the  flow  angle  at  the  trailing  edge.  The  apparatus  and  its 
installation  is  shown  in  figure  3. 


Experimental  Results 

The  experimental  results  are  presented  in  terms 
of  mean  profiles  for  the  velocities  and  Reynolds 
stresses,  and  Fourier  amplitudes  and  phase  lags  for 
these  quantities.  Only  fundamental  harmonic  (twice 
the  elliptic  cylinder  rotation  rate)  data  is 
presented. 

Figure  9  shows  two  cases  of  steady  flow  mean 
velocity  profiles  at  x/c  =  .94,  Ot  =  0,  correspond  inn 
to  horizontal  and  vertical  elliptic  cylinder  orien¬ 
tations.  Figure  10  shows  the  identical  quantities 
for  ct  =  10°.  The  difference  between  the  velocities 
at  the  horizontal  and  vertical  orientations  is  the 
quasi-steady  amplitude  of  the  system.  Note  that 
the  profiles  with  a  steady  vertical  elliptic  cylin¬ 
der  are  closer  to  separation,  and  correspond  to 
much  thicker  boundary  layers. 

The  graph  of  mean  u  velocity  versus  distance  is 
shown  irv  Fig.  11  for  the  case  where  the  elliptic 
cylinder  was  rotating  at  k  =  0.5.  Three  curves  are 
plotted,  corresponding  to  three  local  mean  pressure 
gradients,  f  -  “.0025,  -.014,  -.041.  The  same  geo¬ 
metric  situation  is  replotted  in  Fig.  12  for  k  =  6. a. 
The  increasing  pressure  gradient  makes  the  profile 
fuller  near  the  wall  and  more  concave  near  the  outer 
edqe  of  the  boundary  layer. 

The  mean  velocity  profiles  at  approximately  the 
same  local  mean  pressure  gradient  of  T  ^  -.014, 

for  4  reduced  frequencies  k  =  0.5,  1.0,  2.0  and  6.4 
are  shown  in  Fig.  13,  Except  for  the  k  =  6.4  case 
for  which  T  is  slightly  smaller,  the  differences 
between  the  profiles  are  no  greater  than  the  experi¬ 
mental  error.  This  slight  difference  in  Buri's 
pressure  gradient  parameter  is  probably  due  to  the 
increased  mean  circulation  resulting  from  higher 
cylinder  rotation  rate.  The  effect  of  k  on  the 
mean  pressure  distribution  when  cx  -  0  was  shown  in 
Fig.  4. 

However,  for  the  higher  pressure  gradient  case, 

T  =  -.030,  seen  in  Fig.  14,  differences  due  to  re¬ 
duced  frequency  are  clear.  Increasing  reduced 
frequency  has  a  similar  effect  on  the  profiles  as 
reduced  pressure  gradient.  The  inflection  of  the 
>rofiles  is  reduced,  thus  the  profile  looks  less 
Like  a  separation  profile. 

Ensemble  averaged  velocity  profiles  are  plotted 
in  Fig.  15  for  mean  pressure  gradient  T  =  -.014  and 
reduced  frequency  0.5.  The  profiles  give  the  en¬ 
semble  averaged  velocity  versus  distance  for  several 
representative  times  in  the  rotational  cycle.  For 
this  pressure  gradient  all  profiles  are  without 
inflection  points  at  all  times.  The  reduced  unsteady 
amplitudes  at  high  frequencies  are  seen  in  Fig.  16, 
for  T  =  -.010  and  k  =  6.4.  All  profiles  remain 
well  behaved. 

Fig.  17,  for  k  =  0.5  and  V  =  -.041  shows  quite 
different  character,  as  the  profile  shifts  from 
large  inflection  and  incipient  separation  to  no  in¬ 
flection  during  the  rotation  cycle.  Figure  1h  illus¬ 
trates  higher  frequency  k  -  6.4  at  V  =  --041.  Note 
a  moderately  inflected  profile  is  maintained  through¬ 
out  the  cycle  probably  due  to  smaller  unsteady 
amplitudes.  These  difference',  with  time  f-r  the 
low  frequency  case  are  quite  similar  to  those  ex¬ 
hibited  by  the  steady  profiles  (Figs.  9  and  10). 


The  unsteady  velocity  data  of  last  four  figures 
may  also  be  reduced  to  distributions  of  Fourier 
amplitude  and  phase  lag  across  the  boundary  layer. 
Figure  19  presents  the  amplitude  of  the  fundamental 
harmonic,  normalized  by  the  amplitude  of  the  same 
harmonic  of  the  external  flow,  for  the  same  cases 
as  Fig.  13,  T  ^  -.014,  k  =  0.5,  2.0,  6.4.  The  high 
scatter  seen  for  k  —  6.4  is  a  result  of  an  insuffic¬ 
ient  number  of  averages  for  the  relatively  small 
velocity  amplitude  present.  The  most  striking 
feature  of  these  plots  is  the  large  value  of  the  un¬ 
steady  amplitude  overshoot,  which  ranges  from  2.4 
at  k  *  0.5  to  1.7  at  k  «  6.4.  As  will  be  discussed 
below,  these  values  are  quite  similar  to  previous 
results  for  similar  mean  profiles.  The  height  at 
which  this  maximum  occurs  also  decreases  with  in¬ 
creasing  reduced  frequency. 

Figure  20  shows  the  same  quantities  at  a  higher 
pressure  gradient,  T  =  -.030^  In  addition,  the 
vertical  velocity  amplitude  v  ^.s  also  shown.  As  in 
the  lower  f  case,  the  maximum  u  amplitude  overshoot 
occurs  at  k  =  0.5  with  lower  values  for  k  =  1.0,  2.0 
and  6.4.  However  the  height  of  this  maximum  decreases 
less  rapidly  with  increasing  reduce^  frequency  than 
in  the  data  shown  in  Fig.  19.  The  v  amplitudes  all 
drop  smoothly  to  0  at  the  wall. 

Comparing  Figs.  19  and  20,  the  effect  of  pres¬ 
sure  gradient  is  much  more  apparent  for  lower 
frequencies,  as  the  maximum  overshoot  increases 
from  2.4  times  the  external  amplitude  to  5.6  times 
the  external  amplitude  at  k  =  .5,  while  remaining 
nearly  constant  for  k  =  6.4.  The  height  of  Liiu 
maximum,  at  constant  k,  also  increases  for  higher 
pressure  gradients. 

The  next  property  of  the  velocity  profiles  to 
be  illustrated  are  the  mean  nond intension*  1  Reynolds 
stresses.  Figure  21  shows  u  *uT  for  the  steady  con¬ 
ditions  of  Fig.  9.  Note  the  large  difference  in 
profile  shape  between  the  high  and  low  pressure 
gradient  states.  Increased  pressure  gradient  in¬ 
creases  and  broadens  the  maximum  u 1 u r ,  removing 
the  sharp  peak  near  the  wall.  The  time  average 
value  of  uTur  in  the  presence  of  the  unsteady  flow 
is  shown  in  Fig.  22  for  T  =  -.014  and  V  -  -.010. 
Because  this  data  was  recorded  in  analog  forir ,  the 
experimental  scatter  was  larger  than  for  data  re¬ 
corded  digitally,  such  as  in  Fig.  24.  Never- the- 
less,  these  profiles  do  show  the  peak  very  close  to 
the  surface,  such  as  was  observed  for  Patel's  zero 
pressure  gradient  results^,  and  in  the  low  T  steady 
data  of  Fig.  21.  Values  of  u’u'  to  compare  with  the 
higher  pressure  gradient  steady  case  of  Fig.  10  are 
shown  in  Fig.  23.  The  position  of  the  maximum  has 
moved  further  from  the  wall,  and  the  maximum  has  in¬ 
creased  to  .025  for  the  nearly  .separated  case.  Data 
taker,  with  the  cross  wire  and  the  digital  data  ac- 
(juisit  ion  system,  which  include  measurement  of  all 
three  fjuai.t it ios ,  u’u*,  v’v*,  and  u'v',  are  shown  in 
Fig.  24.  For  this  higher  j^ressu^e  gradient 
(!’  ~  -.030)  case ,  both  u'u'  and  v'v'  appear  indepen¬ 
dent  of  k,  win  If*  u'v'  is  smaller  for  k  -  0.5  than 
for  higher  frequencies.  The  major  effect  of  the 
increased  i  re q sure  gradient  is  to  remove  the  near 
wall  peak  in  u'u',  and  to  increase  the  maximum  value 
by  a  factor  of  three,  from  .005;:  to  .015.  This 
offv.g  j«;  ; . ; ,  1 1  a  r  to  the  effect  of  the  pressure 

,r.,<ii»i,t  «>n  ti;«  ur.stf.idy  amplitude-  discussed  above. 


Finally,  the  amplitude  of  the  fundamental  har¬ 
monic  of  <u’u'>  for  T  =  -.030  and  k  =  0.5,  2.0  and 
6.4  is  shown  in  Fig.  25.  Two  peaks  are  seen,  the 
largest  somewhat  above  the  position  of  maximum  mean 
Reynolds  stress  u‘u 1 ,  and  the  second  close  to  the 
wall.  All  amplitudes  decrease  sharply  as  frequency 
is  increased,  for  example,  the  maximum  of  <u'u'>  is 
.0085  at  k  =  .5,  but  only  .0003  at  k  =  6.4.  This  is 
probably  due  to  both  the  decreased  external  <u> 
amplitude  and  the  decreased  maximum  ratio  of  <u>  to 
the  external  value.  Figure  26  shows  amplitudes  for 
the  three  quantities  <u'u'>,  <u'v‘>,  and  <v'v'>  for 
the  case  k  =  0.5,  T  =  -.030.  The  characteristic 
double  maximum  is  seen  in  each  curve,  however,  the 
height  of  the  upper  maximum  differs. 


Discussion 

As  mentioned  previously,  for  high  reduced  fre¬ 
quencies,  k  >  2.0,  the  unsteady  u  profile  tends  to 
become  independent  of  pressure  gradient,  and  only 
weakly  dependent  on  frequency.  Figure  27  shows  L':.e 
ratio  of  the  maximum  unsteady  amplitude  to  the  ex¬ 
ternal  amplitude  for  all  of  the  boundary  layers 
presently  studied,  with  -.002<F<  -.043.  While  large 
differences  exist  for  k  ^  2,  the  higher  k  values 
definitely  approach  a  common  value  of  approximately 
1.8. 

Figure  27  also  shows  results  from  previous  high 
frequency  data.  In  computing  a  reduced  frequency 
for  flat  plate  experiments,  a  length  equal  to  1/2 
the  distance  from  the  leading  edge  or  from  the 
boundary  layer  trip  is  used.  These  data,  although 
generally  of  somewhat  lower  reduced  frequency,  do 
support  the  relative  independence  of  the  unsteady 
profile  on  pressure  gradient  and  frequency.  This 
result  parallels  theoretical  predictions  of 
Liqhthill  and  Lin  for  laminar  boundary  layers  '  . 

While  not  presented  here  in  detail,  phase  dis¬ 
tributions  of  u  across  the  boundary  layer  were 
generally  small  and  in  agreement  with  those  of 
Cousteix7.  That  is  over  the  outer  part  of  the  boun¬ 
dary  layer  phase  lags  with  respect  to  the  outer  flow 
of  10  to  20  degrees  were  observed,  with  the  position 
of  maximum  lag  corresponding  to  that  of  maximum 
amplitude.  As  the  airfoil  surface  was  approached 
the  phase  changed  to  a  small  lead  (5-10°)  over  the 
external  flow.  Due  to  the  relative  sizes  of  the 
present  hot  wire  nrobe  and  boundary  layer ,  data  in 
the  Viscous  sublayer  was  not  obtained. 

The  ensemble  averaged  profiles  were  normally 
the  only  data  recorded.  However,  for  selected  data 
points,  15-25  cycles  of  raw  data  were  stored.  One 
result  of  the  analysis  of  this  data  was  an  estimate 

of  the  magnitude  of  the  coupling  terms  u'u,  v'u, 

%  n,  -i 

u'v,  v'v.  These  terms  are  commonly  neglected  . 

Based  on  the  present  data  the  coupling  terras  were 
found  to  be  at  least  one  order  of  magnitude  less 
than  the  equivalent  turbulent  terms,  u'u',  u'v', 
v'v* .  It  is  likely  that  if  more  cycles  were 
analyzed  for  this  purpose  the  ratio  would  decrease 
even  further.  From  a  theoretical  viewpoint  such 
coupling  terms  should  vanish  if  the  ensemble  average 
is  taken  at  the  same  point  in  each  cycle. 


The  unsteady  portion  of  the  ensemble  averaged 
Reynolds  stresses  as  shown  in  Fig.  25  as  a  Fourier 
amplitude  is  a  different  way  of  illustrating  the 
phenomena  observed  previously  in  Reference  7. 
Cousteix,  et  al  ,  show  a  time  dependent  variation 
in  the  maximum  <u’v'>.  It  varies  from  -.0055  to 
-.002,  well  within  the  range  of  the  present  data, 
i.e.  <u'v'>  t  -  .00375  +  .00175  cos  from  the 
Fourier  analysis  cited  above. 

Figure  28  presents  another  significant  feature 
of  the  ensemble  averaged  Reynolds  stress:  the  phase 
lag  across  the  boundary  layer  for  a  typical  low 
frequency,  high  pressure  gradient  (and  therefore  high 
amplitude)  case.  The  phase  of  <u'u'5  is  seen  to  be 
about  180°  different  from  that  of  <u>  over  the  outer 
boundary  layer,  with  a  shift  taking  place  at  approxi¬ 
mately  one  mean  momentum  thickness  from  the  wall,  to 
becoming  within  10°  of  being  in  phase  with  <u>.  In 
terms  of  the  mean  profile,  <u>  and  <u'u’>  are  out  of 

phase  when  and  ^  (u'u')  are  of  opposite  sign, 

and  in  phase  when  they  are  of  the  same  sign. 


CONCLUSIONS 

A  series  of  experiments  have  been  performed  to 
determine  the  characteristics  of  the  turbulent  boun¬ 
dary  layer  of  a  NACA  0012  two  dimensional  airfoil 
subject  to  an  oscillating  external  flow.  Ensemble 
averaged  velocities  and  Reynolds  stresses  were 
measured  over  a  wide  range  of  mean  pressure  grad¬ 
ients  and  reduced  frequencies,  from  nearly  quasi 
steady  to  signif icantly  greater  than  a  reduced 
frequency  of  unity,  and  from  a  mildly  adverse 
gradient  to  one  producing  incipient  separation. 

Mean  profiles  were  found  to  be  quite  indepen¬ 
dent  of  reduced  frequency  for  moderate  adverse 
pressure  gradients,  but  for  steeper  adverse  pres¬ 
sure  gradients,  the  mean  profile  was  less  like  a 
separation  profile  as  the  frequency  was  increased. 

Periodic  velocity  profiles  had  large  increases 
in  amplitude  in  the  boundary  layer  as  compared  to 
external  flow,  but  small  phase  differences  from 
that  flow.  For  low  frequencies,  the  profiles  were 
highly  dependent  on  the  mean  pressure  gradient, 
while  the  high  frequency  profiles  appeared  to 
approach  a  universal  profile  independent  of  both 
frequency  and  pressure  gradient. 

Ensemble  averaged  Reynold  stresses  had  a  much 
stronger  dependence  on  the  mean  pressure  gradient 
than  on  reduced  frequency,  with  both  high  and  low 
pressure  gradient  data  being  similar  to  previous 
results  for  roughly  similar  conditions. 

Correlations  between  the  periodic  and  non¬ 
periodic  velocity  components  (u'u)  were  found  to  be 
at  least  an  order  of  magnitude  smaller  than  the 
mean  nonperiodic  values  (u’u' ) .  Theoretically  this 
value  should  be  zero.  The  common  neglect  of  these 
quantities  appears  to  be  justified. 

A  Fourier  component  of  the  ensemble  averaged 
Reynolds  stress  was  extracted  from  the  data  and  was 
found  to  depend  on  the  amplitude  of  the  periodic 
velocities,  and  on  the  relative  sign  between  the 
normal  gradients  of  the  mean  velocity  and  mean 
Reynolds  stress. 
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Fig .  1 


External  Velocity  Distribution  on 
Flat  Plate 


Fig.  2  Mean  Velocity  Profile  on  Flat  Plate 


Fig.  3  Schematic  Diagram  for  Unsteady 
Airfoil  Test 

1)  Test  Section,  2)  Airfoil, 

3)  Rotating  Elliptic  Cylinder, 

4)  Drive  Motor,  5)  Sidewalls, 
6)  Pitot-Static  Probe 


Fig.  4  Mean  Airfoil  Pressure  Distributions, 
a  =  o° 


Fig.  S  Mean  Airfoil  Pressure  Distributions, 

a  «  10° 
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Fig.  6  Fundamental  Harmonic  Pressure  Amplitude 
Distributions 
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Fig.  7  Fundamental  Harmonic  Pressure  Phase 
uag  Distributions 
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Fig.  8  Instrumentation  Block  Diagram 
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Fiq.  0  Mean  Velocity  Profiles  for  Steady 
Elliptic  Cylinder,  a  -  0° 
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Fig.  10  Mean  Velocity  Profiles  for  Steady 
Elliptic  Cylinder,  a  =  10° 


Fig.  11  Variation  in  Mean  Velocity  Profiles  Due 
to  Mean  Pressure  Gradient  at  Constant 
Reduced  Frequency,  k  =  0.5 


Fiq.  12  Variation  in  Mean  Velocity  rrofiles  Due 
to  Mean  Pressure  Gradient  at  Constant 
Reduced  Frequency  k  =  6.4 
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Fig.  21  Mean  Reynolds  Stress,  u'u'.  Steady 
Elliptic  Cylinder  and  a  =  0° 


Fig.  24  Mean  Reynolds  Stresses,  u'u1,  v'v1  and 
u’v'.  Rotating  Elliptic  Cylinder  and 
a  =  10° 
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